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Recent work on many particle system reveals the existence of regular collective perturbations 
corresponding to the smallest positive Lyapunov exponents (LEs), called hydrodynamic Lyapunov 
modes. Until now, however, these modes are only found for hard core systems. Here we report new 
results on Lyapunov spectra and Lyapunov vectors (LVs) for Lennard-Jones fluids. By considering 
the Fourier transform of the coordinate fluctuation density u^ a \x,t), it is found that the LVs with 
A ~ are highly dominated by a few components with low wave-numbers. These numerical results 
provide strong evidence that hydrodynamic Lyapunov modes do exist in soft-potential systems, 
although the collective Lyapunov modes are more vague than in hard-core systems. In studying the 
density and temperature dependence of these modes, it is found that, when the value of Lyapunov 
exponent A^"' is plotted as function of the dominant wave number k max of the corresponding LV, all 
data from simulations with different densities and temperatures collapse onto a single curve. This 
shows that the dispersion relation vs. k ma x for hydrodynamical Lyapunov modes appears to be 
universal irrespective of the particle density and temperature of the system. Despite the wave-like 
character of the LVs, no step-like structure exists in the Lyapunov spectrum of the systems studied 
here, in contrast to the hard-core case. Further numerical simulations show that the finite-time LEs 
fluctuate strongly. We have also investigated localization features of LVs and propose a new length 
scale to characterize the Hamiltonian spatio-temporal chaotic states. 

PACS numbers: 05.45.-a,05.20.-y,02.70.Ns,05.20. Jn,05.45. Jj 



I. INTRODUCTION 

One of the most successful theories in modern science 
is statistical mechanics, which allows one to understand 
the macroscopic (thermodynamic) properties of matter 
from a statistical analysis of the microscopic (mechani- 
cal) behavior of the constituent particles. In spite of this, 
using certain probabilistic assumptions such as Boltz- 
mann's Stosszahlansatz renders the lack of a firm foun- 
dation of this theory, especially for non-equilibrium sta- 
tistical mechanics. Fortunately, the concept of chaotic 
dynamics developed in the 20th century Q is a good 
candidate for complying with these difficulties. Instead 
of the probabilistic assumptions, the dynamical instabil- 
ity of trajectories can provide the necessary fast loss of 
time correlations, ergodicity, mixing and other dynamical 
randomness 2}. It is generally expected that dynamical 
instability is at the origin of macroscopic transport phe- 
nomena and that one can find certain connection between 
them. In the past decade, some beautiful theories in 
this direction were already developed. Examples are the 
escape-rate formalism by Gaspard and Nicolis 0, Q and 
the Gaussian thermostat method due to Nose, Hoover, 
Evans, Morriss and others 0,0,0' where the Lyapunov 
exponents were related to certain transport coefficients. 

Very recently, molecular dynamics simulations on 
hard-core systems revealed the existence of regular col- 
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lective perturbations corresponding to the smallest pos- 
itive Lyapunov exponents (LEs), named hydrodynamic 
Lyapunov modes This opens a possible new way 
for the connection between Lyapunov vectors, a quan- 
tity characterizing dynamical instability of trajectories, 
and macroscopic trans por t properties. A lot of work 
UdmElliQElIl has been done to identify this 
phenomenon and to find out its origin. It is commonly 
thought that the appearance of these modes is due to the 
conservation of certain quantities in the systems studied 
flfl ITU Il2l Hal IT3 | . A natural consequence of this expec- 
tation is that the appearance of such modes would not 
be an exclusive feature of hard-core systems and should 
be generic to a large class of Hamiltonian systems. How- 
ever, till now, these modes were only identified in the 
computer simulations of hard-core systems 0, 0, 0] . 

In this work, we report on new results about a Id sys- 
tem with Lennard-Jones interaction. Although the iden- 
tification of regular hydrodynamic Lyapunov modes by 
naked eye is difficult for soft potential systems our 
new technique based on a spectral analysis of LVs shows 
strong evidence that hydrodynamic Lyapunov modes do 
exist in this case. The influence of density and tempera- 
ture changes is studied in detail. The dispersion relation 
for hydrodynamic Lyapunov modes, the dominant wave- 
number as function of the corresponding LEs is found 
being universal for all densities and temperatures. 

Furthermore we study the localization properties of 
LVs. Based on the extensive nature of LVs with A « 0, 
we propose a new length scale to characterize a spatio- 
temporal chaotic Hamiltonian system. It is expected 
that this new quantity will be useful for the task of dis- 
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tinguishing different spatio-temporal chaotic states and 
characterizing transitions among them. This is an im- 
portant open question in the study of spatial-temporal 
chaos IH. 



II. MODEL 

In this study we use a ID Lennard-Jones system with 
Hamiltonian 
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The interaction potential among particles is of the form: 
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FIG. 1: Time evolution of temperature T = (mv 2 ) and total 
energy. The nearly constant state variables show that the 
system has reached already a stationary state. The parameter 
setting used here is: N = 100, L = 1000 and T = 0.2. 
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with V c = 4e [(a/r c ) 12 - (a/r c ) e ] . The potential is trun- 
cated in order to lower the computational burden. Other 
types of potential with smoothed force at the truncation 
point were also simulated to check the influence on the 
results given below. No qualitative difference is found 
between them. 

The system is integrated using the velocity form of 
the Verlet algorithm with periodic boundary conditions 
|18|. In our simulations, we set m = l,<r = l,e=l 
and r c = 2.5. All results are given in reduced units, i.e., 
length in units of a, energy in units of e and time in 
units of (mcr 2 /48e) 1 / 2 . The time step used in the molec- 
ular dynamics simulation is h — 0.008. The standard 
Gram-Schmidt re-orthonormalization algorithm 0, [2(j 
is used to calculate the local dynamical instability of 
the systems studied. The time interval for periodic re- 
orthonormalization is 30ft. to 100ft. Throughout this pa- 
per, the particle number typically is denoted by N, the 
length of the system by L and the temperature by T. 



III. NUMERICAL RESULTS 

A. The stationary state 

In the numerical calculation of the Lyapunov instabil- 
ity of a many-body system [2l| , there are some important 
time scales to be kept in mind: the first one is the time for 
a many-body system to relax to a stationary state, which 
guarantees that quantities measured afterwards are not 
for a transient state; the second is the time for the set 
of Lyapunov vectors to relax to their right orientations 
since offset vectors are usually selected randomly at the 
beginning; the third is the time used to count LEs and 
LVs, which should be long enough to ensure that the tra- 
jectory wanders all over the attractor. For a many-body 
system like the one studied here, these time scales can 
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FIG. 2: a) Particle positions x% vs. index of particles i and 
b) pair distribution function G(r) for the space among all 
particles (upper panel) and nearest neighbors (lower panel) for 
the stationary state shown in Fig. Q (see 24j for the definition 
of G(r)). The sharp peaks in G(r) imply that the state is a 
broken-chain one. 



be extremely long due to the large number of degrees of 
freedom involved plE^. 

The time evolution of state variables like temperature 
T and total energy for a case with parameter setting 
N = 100, L = 1000 and T = 0.2 is shown in Fig. □ 
In the beginning of our molecular dynamics simulation, 
particles are placed randomly in the interval [0, L] . Their 
velocities are chosen randomly from a Boltzmann distri- 
bution. In order to equilibrate the system, it is coupled 
to a stochastic heat bath with given temperature T, i.e., 
every 500 steps the velocities of the particles are replaced 
with velocities that were drawn from a Boltzmann dis- 
tribution corresponding to that temperature. This was 
done for a time period of length t eq , which is longer than 
the relaxation time of the system at this temperature. 
After the equilibration procedure, the system is allowed 
to evolve with constant total energy, i.e., without the 
heat bath, for a time period of the same length as t eq , in 
order to be sure that the system is already in a stationary 
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state at given temperature T. In Fig. ^ the period with 
thermal bath is omitted and only the part of evolution 
with constant total energy is shown. The nearly constant 
value of temperature means that the system has already 
reached a stationary state and one can start the calcu- 
lation of the Lyapunov instability of the system. The 
pair distribution function G(r) shown in Fig. [21 tells us 
that the stationary state for T — 0.2 is a broken-chain 
state with short range order. This is generic for a Id 
Lennard- Jones system with not too high density p5j . 



B. Benettin method using Gram-Schmidt 
orthogonalization 

The standard method invented by Benettin and Shi- 
mada 0, |2(j is the most efficient one to calculate the 
Lyapunov exponents and Lyapunov vectors of large sys- 
tems. Here N x N linearized equations for offset vectors 
in the tangent space were integrated simultaneously with 
the set of N nonlinear equations for the reference trajec- 
tory. Offset vectors were periodically re-orthonormalizcd 
using the Gram-Schmidt algorithm. The resulting rescal- 
ing factors measure the expansion or contraction rate of 
offset vectors in certain directions. Averaging the loga- 
rithm of them for a period r gets what are called finite- 
time Lyapunov exponents A T . The limit A = A T _>+oo 
is what is usually called Lyapunov exponent. The value 
of finite-time LEs A T depends on the trajectory segment 
where it is calculated and usually it fluctuates as the 
segment moves along the trajectory. However A +oc is 
time independent and unique for an ergodic system. In 
this sense, A+oc is a global quantity in characterizing 
the system attractor while the finite-time LEs are lo- 
cal quantities which contain more detailed information 
about the dynamics. The offset vectors just after re- 
orthonormalization are called Lyapunov vectors. It is a 
local quantity in characterizing the system attractor sim- 
ilar to the finite-time LEs. 

Another point to be noted is that Lyapunov vectors 
obtained using Benettin's method are always mutually 
orthogonal while the local unstable and stable directions 
are in general not orthogonal. In this sense, these are two 
different sets of vectors. They are also different from the 
one in the multiplicative ergodic theorem 26] . Lyapunov 
vectors obtained in the standard way can at least repre- 
sent the most unstable direction in a certain subspace 
and they already contain a lot of important information 
about the dynamical instability in tangent space. We will 
rely on them to continue our study in this paper. 



C. Lyapunov exponents with wild fluctuations 

The Lyapunov spectrum for the case with N = 100, 
L = 1000 and T = 0.2 is shown in Fig. El Here only half 
of the spectrum is shown since all LEs of Hamiltonian 
systems come in pairs according to the conjugate-pairing 
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FIG. 3: Lyapunov spectrum for the state shown in Fig. |51 
Enlargement of the part in the regime A^ a ' ~ is shown in 
the inset. It is the result of an average over 58 samples with 
different initial conditions and for each sample the averaging 
period is 4 x 10 e h. Here no step- wise structure exists in con- 
trast to the case of hard-core system. This is the typical result 
for our soft potential system. 




FIG. 4: Distribution of the finite time Lyapunov exponent 
Af where r is equal to the period of re-orthonormalization 
and the index of LEs a is equal to 90, 94 and 98 respec- 
tively. The strong fluctuations of are one of the possible 
reasons for the disappearance of the step-wise structures in 
the Lyapunov spectrum. The parameter setting used here is: 
N = 100, L = 1000 and T = 0.2. 



rule 
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|28| . From the enlargement shown in the inset 
of Fig. |2| for the part near A 1 -"' w 0, one can not see any 
step- wise structure in the Lyapunov spectrum in contrast 
to the case of hard-core systems 10] . This is the typical 
result obtained for our soft potential system. 

The fluctuations in local instabilities is demonstrated 
by the distribution of finite-time LEs. In Fig. 01 
such distributions for some LEs in the regime A w 
are presented. Fluctuations of the finite time Lya- 
punov exponents are quite large compared with the 
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FIG. 5: Normalized Lyapunov exponents A^/A' 1 ' with 
L = 200, 400, 1000 2000, and 4000 respectively. The Lya- 
punov spectrum becomes more and more bended as the parti- 
cle density p — N / L is decreased. This implies the separation 
of two time scales. Here TV = 100 and T — 0.2. 



FIG. 6: u (a) (x,t) for LVs with index a = 1, 10, 90, and 95 
respectively. Notice that the former two are more localized 
while the latter two are more distributed. The parameter 
setting used here is: N = 100, L = 1000 and T = 0.2. 



difference between their mean values, i.e., tr(Ar ) = 

\/(xi a)2 ) - (\ { T a) ) 2 > |A<°> - A(« +1 )|. Here (• • • ) means 
time average. The strong fluctuations in local instabili- 
ties is one of the possible reasons for the disappearance 
of the step-wise structures in the Lyapunov spectra, ft 
could also cause the mixing of nearby Lyapunov vectors. 
The mixing may be at the origin of the intermittency ob- 
served in the time evolution of the spatial Fourier trans- 
form of LVs (see Sect. MI E(l . 



D. Bending of Lyapunov spectrum with decreasing 
particle density 

We studied also the influence of the particle density 
on the Lyapunov spectrum by increasing the length of 
the system with the particle number N kept fixed. As 
can be seen in Fig|SJ the Lyapunov spectrum becomes 
more and more bended with increasing L. For the case 
of L = 4000, the spectrum can already be unambiguously 
divided into two regimes: In the upper regime, Lyapunov 
exponents decrease more quickly with increasing index 
than in the lower regime. This bending of the Lyapunov 
spectrum was related to the separation of two time scales 
in dilute particle systems |2^. We propose that one is 
the time of local collision events, and the other is due to 
the collective motion of the particles. For a system with 
high density, the collisions are so frequent that there are 
strong correlation between consecutive collisions and one 
can no longer separate them from each other. The col- 
lisions themselves contribute to the collective motion of 
the system. Therefore no time scale separation happens 
here and the LEs decrease gradually. 



E. Spatial structure of LVs with A (q) « 
1. Coordinate fluctuation density (CFD) 

Another quantity used to characterize the local insta- 
bility of trajectories are Lyapunov vectors ST^ a \ which 
represent expanding or contracting directions in tangent 
space. In the study of hard-core systems, Posch et al. 
found that coordinate part of the Lyapunov vectors cor- 
responding to A w are of regular wave-like character 
M ^3 • They are referred to as hydrodynamic Lya- 
punov modes. We are searching here for the counterpart 
of these modes in our soft-potential system. 

Remember that each of the LVs is consist of two parts: 
the displacement Sxi in coordinate space and Svi in mo- 
mentum space. In past studies of hydrodynamic Lya- 
punov modes in hard-core systems, only the coordinate 
part 5xi was considered. This is due to an interesting 
feature of LVs found in Ref.^2] which says that the an- 
gles between the coordinate part and the momentum part 
are always small, i.e, the two vectors are nearly paral- 
lel. Therefore, it is already sufficient to use only Sxi 
for studying ST. For our soft potential system, we find 
that the angles between the coordinate part and the mo- 
mentum part are no longer as small as in the hard-core 
systems. However, we will still follow the tradition to 
study the coordinate part of LV first and come to the 
momentum part afterwards. 

Analogous to the definition of microscopic densities in 
hydrodynamics |24|. we define a quantity called coordi- 
nate fluctuation density 

N 

u^(x,t) = J2$4 a) -S(x-Xi). (3) 

i=l 

Profiles of u^ a \x,t) for some typical LVs of the 
Lennard-Jones system are presented in Fig. It can be 
seen that u^(x,t) for LVs corresponding to the largest 
Lyapunov exponents are highly localized, for example 
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FIG. 7: Time evolution of m (95) (x, t). No clear wave structure 
can be detected. The parameter setting used here is: N = 
100, L = 1000 and T = 0.2. 

u^'(x,t) and u^ 10 ^(x,t), while those for LVgo and LV95 
are more distributed. The study on the localization of 
LVi is of long history dESEllMElllllllllland 
it was related to defect events in simulations of Benard 
convection .'!.'( . We leave the discussion on this point to 
Scc lHI Fl The temporal evolution of u^ 95 \x,t) is shown 
in Fig. [7| in order to make the possibly existing wave- 
like structure more evident. A wave structure however, 
cannot unambiguously be detected here with the naked 
eye. 

2. Spatial power spectrum of CFD and intermittency in its 
time evolution 

Now we turn to the spatial Fourier transform of 
u^^Xjt), which reads 

u^\t) = J u^(x,t)exp(~ikx)dx 

N M (4) 
= 5Xj ■ exp[— ik ■ Xj(t)] 

In previous studies, in order to make the wave struc- 
ture more obvious, certain smoothing procedures in time 
or space were applied to the Lyapunov vectors. For a Id 
hard-core system with only a few particles, this proce- 
dure has been shown to be quite useful in identifying the 
existence of hydrodynamic Lyapunov modes 15]. The 
success of this strategy relies on the fact that some of 
the Hydrodynamic Lyapunov modes (transverse modes) 
of hard-core systems are stationary |8( • Therefore, time 
averaging can indeed suppress the noise component and 
make the long wave-length modes more significant. For 
a soft potential system, all the Lyapunov vectors are not 
stationary due to the random mixing among them. Es- 
pecially, for our one dimensional system studied here, no 
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FIG. 8: Intermittent behavior of the peak wave-number k, 
and spectral entropy H s (i) for the spatial Fourier spectrum 
of t/ 95 ^ (x,t). The parameter setting used here is: N = 100, 
L = 1000 and T = 0.2. 

transverse modes but only longitudinal Lyapunov modes 
exist. In consequence of this, the smoothing procedure 
is no longer very helpful for detecting the hidden regular 
modes and can even damages them |lul . Here we apply 
the spatial Fourier transformation to the instantaneous 
quantity u^ aS) (x,t) instead. The algorithm offered espe- 
cially for unevenly distributed data is very suitable for 
our case |35j . Furthermore, we take the long time average 
(and ensemble average) of the spatial Fourier spectrum 

s^(k,t) = \u^(t)\ 2 (5) 

since it is expected that in 5i« (k) = (s£u }(k, t)) the con- 
tribution of stochastic fluctuations will be averaged out 
while the information about the collective modes will be 
left and accumulated. The following results show that 
this technique is quite successful in detecting the hidden 
collective modes. 

The time evolution of the spatial Fourier spectrum 
Suu\k,t) for Lyapunov vector No. 95 is investigated in 
Fig. IS] as an example. Two quantities are recorded with 
time going on. One is the peak wave- number which 
marks the position of the highest peak in the spectrum 
Suu (k,t) (see Fig. |5J. The other is the spectral entropy 
H s (t) 36], which measures the distribution property of 

the spectrum s««(fc,t). It is defined as: 

H s (t) = -J2 *) ^ { uJ(h,t) (6) 

ki 

A smaller value of H s (t) means that the spectrum 

Sum (k, t) is highly concentrated on a few values of k, 
i.e., these components dominate the behavior of the LV. 
Both of these quantities behave intermittently as shown 
in Fig. |S| Large intervals of nearly constant low values 
(off state) are interrupted by short period of bursts (on 
state) where they experience large values. Details of typ- 
ical on and off states are shown in Fig. [3] It can be seen 
that the off state is dominated by the low wave-number 
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FIG. 9: a) Variation of the peak wave number k* with time. 
b),c) Two typical snapshots of LV95, off and on state at t — 44 
and 176 respectively. d),e) their spatial Fourier transform. 
The spectrum for the off state has a sharp peak at small 
fc* while that for the on state has no dominant peak. The 
parameter setting used here is: TV = 100, L = 1000 and 
T = 0.2. 



components (see the sharp peak at low wave-number 
while the on state is more noisy and there are no signif- 
icant dominant components. This intermittency in the 
time evolution of the spatial Fourier spectrum of LVs is 
a typical feature of soft potential systems. It is conjec- 
tured that this is a consequence of the mixing of nearby 
LVs caused by the wild fluctuations of local instabilities. 
One time scale can be extracted from the intermittency, 
i.e., the mean duration r a ff of the off state. We con- 
jectured that this time scale is related to the life time 
of hydrodynamics Lyapunov modes. If the off state was 
viewed as a pure mode, the time r Q ff is just the aver- 
age life-time of such modes. Due to the mutual interac- 
tion among modes, the hydrodynamics Lyapunov modes 
in the soft potential systems are only of finite life-time. 
Further numerical work is needed in this direction to get 
more direct evidence. 

In Fig llOl the time averaged spectral entropy (H s (t)) 
is plotted against the index of the LVs. It increases grad- 
ually as the index decreases from N — 2. This means 
that LVs corresponding to smaller positive LEs are more 
localized in Fourier space, i.e., they have more wave- like 
character than those corresponding to larger LEs. 



3. Dispersion relation of hydrodynamical Lyapunov modes 

Now we consider the time averaged spatial Fourier 
spectrum S^H(k) of LVs. Two cases with L = 1000 and 
2000 are shown in Fig^] It is not hard to recognize the 
sharp peak at A ps in the contour plot of the spectrum. 
In increasing the Lyapunov exponents, the peak shifts to 
the larger wave number side. A dashed line is plotted 
to guide eyes how the wave number of the peak k max 
changes with \( a \ To further demonstrate this point, 
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FIG. 10: Time-averaged spectral entropy (H s (t)} vs. index of 
LVs. The gradual decrease of (H s (t)) from a m means that 
LVs corresponding to smaller positive LEs are more localized 
in Fourier space, i.e. they have more wave- like character those 
that corresponding to larger LEs. The parameter setting used 
here is: TV = 100, L = 1000 and T = 0.2. 



in Fig ll2l the value of the Lyapunov exponent A^ Q ^ is 
plotted versus k max of corresponding LVs. We call this 
the dispersion relation of the hydrodynamical Lyapunov 
modes. The numerical fitting of the data shows that for 
A ~ 0, A*"' ~ k max with the exponent 7 ps 1.2. However, 
a linear dispersion with quadratic corrections can not be 
excluded. 

In order to show that the peak in S^u (k) is not a result 
of the highly regular packing of particles in the broken- 
chain state, the static structure function [24| 



S(fc) = / G(r) exp(-ikr)dr 



(7) 



for the case L — 2000 is plotted in the same figure 
as S^}(k), where G(r) is the pair correlation function 
shown in Fig. [3 It can be seen that S(k) is nearly 
constant in the regime k ps 0, the place where a sharp 
peak was observed in Stm(k). The regular packing of 
particles causes the formation of a peak at fc/27r ps 0.9 
in S(k). This corresponds to a tiny peak at the same 
&;-value in S^(k) for those LVs with A ps 0. These facts 
show clearly that the collective modes observed in LVs 
are not due to the regular packing of particles. 

To further demonstrate the properties of these modes, 
in Fig^l k max is plotted versus index of LVs. As can be 
seen, for LVs with A s» 0, i.e. with a pa N, the value of 
kmax is quite small (see the enlargement in Fig. 114(1 . For 
example, for a = 96, 97, and 98, k max = 2tt/L, which is 
the smallest wave number allowed by the periodic bound- 
ary conditions used. Another point to be noticed is the 
step structure in plotting k max as function of a. This is 
similar to the degeneration of wave-numbers found in the 
hard-core case, although the steps here are not so regular. 
In the middle panel of Fig^D the height Su$(k max ) of 
the highest peak in the time averaged spatial spectrum is 
plotted with index of LVs. Except LV No. 99 and No. 100 
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FIG. 11: (Color online) Contour plot of the spectrum 
for L = 1000 and 2000. There is a sharp peak at k m and 
A w 0. To guide eyes, a dashed line is plotted to show how 
the peak wave-number k max changes with the variation of A. 
The sudden jump in k max is marked with an arrow. In total 
58 samples for the case of L = 1000 (10 for L = 2000) are 
used averaging each for a period of 4 x 10 6 h. Here T — 0.2 
and N = 100. 
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FIG. 13: (Upper panel) The wave-number k ma x of the high- 
est peak in the time averaged spatial Fourier spectrum of LVs. 
(Middle panel) The height 5*i"' (k max ) of the highest peak in 
the time averaged spectrum. (Lower panel) The spectral en- 
tropy Hs for the averaged spectrum SiZ > (k). The sudden 
jump in kmax implies the separation of time scales. The pa- 
rameter setting used here is: N = 100, L = 1000 and T = 0.2. 
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FIG. 14: Enlargement of Fig |13l for the part in the regime 
a ~ N. 




FIG. 12: The Lyapunov exponent A (a) is plotted as function 
of the wave-number k ma x of the highest peak in the time 
averaged spatial Fourier spectrum of LVs. The dashed line is 
of the form A' a ' ~ k m 2 ax . The parameter setting used here is: 
N = 100, L = 1000 and T = 0.2. 



for the conserved quantities, Sm (k max ) decreases grad- 
ually in decreasing the index from N — 2. Similar to 
the definition of the spectral entropy for the instanta- 
neous spectrum s^}(k,t) in Eq|^l one can also define a 
spectral entropy Hs for the averaged spectrum S^(k). 
The spectral entropy Hs for S^}(k) is presented in the 
lower panel of Fig ll3l It possesses a minimum at a — 98 
where the Lyapunov exponent is the smallest positive 
one. According to the definition of the spectral entropy, 
the minimum means that the spectrum of LV No. 98 is 
most significantly dominated by a few components. 

All of our results shown above give strong evidence 
that the Lyapunov vectors corresponding to the smallest 
positive LEs in our Id Lennard- Jones system are highly 
dominated by a few components with small wave num- 
bers, i.e, they are similar to the Hydrodynamic Lyapunov 
modes found in hard-core systems. The wave-like char- 
acter becomes weaker and weaker as the value of LE is 
increased gradually from zero. 
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FIG. 15: The Lyapunov spectrum A'"' (upper panel), 
d 2 \/da 2 (middle panel) and k max vs. the index a. The pa- 
rameter setting used here is: N = 100, L = 1000 and T = 0.2. 




FIG. 16: Same as top panel of Fig. |H| but with different 
density (a) p = 1/2, (b) 1/4 and (c) 1/8 respectively. Here 
T = 0.2 and the particle number TV = 100. 
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FIG. 17: (Color online) (upper panel) The wave-number kmax 
of the highest peak in the time-averaged spatial Fourier spec- 
trum of LVs as function of LV index a. (lower panel) the 
Lyapunov exponent A^ a ' as function of kmax- Note that in 
the lower panel, all data from simulations with different den- 
sities and temperatures collapse to a single curve. Fitting of 
the low wave-number part to a power-law function X a ~ k max 
gives 7 » 1.2 ± 0.1. Here N = 100 and T = 0.2. 



4- Separation of time scales 

Another interesting point in Fig |13l is the sudden jump 
in k max at a « 26 which divides the whole set of LVs 
into two groups. It is believed that this sudden jump 
is related with the bending of the LE spectrum and the 
separation of time scales in a dilute system. As shown 
in Fig. 1151 the sudden jump is in the regime where the 
LEs spectrum is most strongly bended, although it is 
not at the exact place where d 2 X^ /da 2 experiences the 
maximal value. Further work is needed to reveal the 
underlying connection between these phenomena. 

5. Influence of density and temperature 

To study how the change in density influences the be- 
havior of LVs, we increase the length L of the system 
from 200 to 4000 with the particle number N kept fixed 
at 100. From the time evolution of fc* shown in Fig ll6l 
one sees that, in increasing the density p = N/L, the 
occurrence of the on-state becomes more frequent, i.e., 



the domination of low wave-number components is much 
weaker. The spatial Fourier spectrum for LVs with LEs 
in the regime A' Q ) ~ 0, however, are always dominated 
by certain low wave-number components irrespective of 
the density (see FigfTTjl. 

An important point is the collapse of data of disper- 
sion relations from simulations with various densities and 
temperatures to a single curve (see lower panel of Fig |T7) l. 
This means, for hydrodynamic Lyapunov modes in our 
system, there is a universal dispersion function A Q (fc) ir- 
respective to the particle density and the system tem- 
perature. Fitting of the data to a power-law function 
X a ~ fc^ oa . says the value of the exponent 7 is 1.2 ± 0.1 
although a linear dispersion with quadratic corrections 
cannot be excluded. 

Another feature of FigEI( u PP er panel) is that the sud- 
den jump in k max disappears as the system density N/L 
is increased. This is consistent with our above discussion 
that the separation of time scales is only significant in 
dilute systems. 
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FIG. 18: Time-averaged participation number vs. index 
of LVs (upper panel) and X^ (lower panel). The parameter 
setting used here is: N = 100, L = 110 and T = 1.6. 
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FIG. 19: Time-averaged participation number p"' and 
pi N ~ 2 ) vs. particle number N. The parameter setting used 
here is: L/N = 1.1 and T = 1.6. 



6. Dynamics of the momentum part 

Now we turn to investigations on the spatial Fourier 
spectrum of the momentum part of LVs. Unfortunately, 
all the spectra are more or less homogeneously dis- 
tributed on all wave- numbers. For all the cases checked, 
no wave-like structure can be identified as for the coordi- 
nate part. One may wander why no mode-like collective 
motion is observed in the momentum part. There are two 
possibilities, one is that the momentum part does contain 
information similar to the coordinate part but it is too 
weak to be detected here due to the strong noise. The 
other is that there is no similarity between the two parts 
at all and regular long wave-length modes exist only in 
the coordinate part. Further work is needed to clarify 
which one is correct. 




12 0.4 0.6 0.8 1 1.2 1.4 1.6 

T 



FIG. 20: a) Time-averaged participation number p^ and 
p( N ~ 2 > vs. temperature T. The parameter setting used here 
is: N = 200 and L = 220. 



F. Localization properties of LVs 

To study the localization of LVs, we employ the par- 
ticipation number, which is defined as: 



r N 



2\2 



L i=l 



(8) 



for a Lyapunov vector (Sx^ , Sv^ ) j37j . This is a stan- 
dard quantity used in the study of disorder-induced local- 
ization , which roughly measures the number of parti- 
cles which contribute to the Lyapunov vector. For the ho- 
mogeneous Lyapunov vector LVn with w 2 = 5xf + 5vf = 
1/N, which corresponds to one of the zero- value LEs, p 
gets its maximal value N . In decreasing the index of LVs, 
LEs become larger and larger. Accompanying to this, 
the participation number p^ decreases as presented in 
Fig. [^J where the variation of the time-averaged value 
of p( a ) versus index a and X^ are plotted. The decrease 
of p( a ) implies that LVs become more and more localized 
with decreasing a. 



One should note that p( N 2 - ) for LVn-2 correspond- 
ing to the smallest positive Lyapunov exponent is signif- 
icantly different from = N for LV N with X^ N ~> = 0. 



In Fig. HI p 



(JV-2) 



40 while p 



(A ; ) 



100. In the 



study of space-time chaos, a commonly used measure for 
quantifying spatio-temporal disorder is the fractal dimen- 
sion D. A spatio-temporal chaotic system is generally 
extensive, with the fractal dimension D in proportion 
to its volume V. According to this, a bounded inten- 
sive quantity named dimension correlation length is de- 
fined £d = limv^oo(D/V)~ 1 / d . Based on the intuitive 
thought that a spatio-temporal chaotic system is com- 
posed of many subsystems and that these subsystem are 
uncorrelated if they are far apart, it is expected that 
£d is proportional to the two-point correlation length, 
which measures the spatial disorder of the system. A par- 
ticular fractal dimension, the Lyapunov dimension Dc 
can be easily obtained using the Kaplan- Yorke formal- 
ism which relates the Lyapunov dimension to Lyapunov 
exponents of the system [39j . For a Hamiltonian system 
like the one studied here, Lyapunov exponents are paired, 
i.e. A'"* 1 + ^( 2dAf - a ) = due to the sympletic structure 
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of the system [271 US • According to the Kaplan- Yorke 
formalism, the Lyapunov dimension for a d-dimensional 
Hamiltonian spatio-temporal chaotic will be Dc = 2dN 
in spite of the details of dynamics. Although Dc de- 
fined in such way is proportional to the volume of the 
system i.e., reflects the extensiveness of the system, its 
value is always a constant irrespective of the temperature 
change. However the two-point correlation length does 
change with the temperature and it even becomes diver- 
gent as a phase transition is encountered. In this sense, 
the Lyapunov dimension Dc and the dimension corre- 
lation length £d are only trivially defined here and not 
good quantities to characterize spatio-temporal Hamilto- 
nian chaos. Here we propose a new length scale based on 
the participation number of LVs: 

. / p (o) \V d 

- ( V • v ) ^ 

where d is the dimension of the physical space and V 
is the system volume which is simply L for our d = 1 
case here. According to the fact we mentioned above, 
£p N 2 ' for the smallest positive Lyapunov exponent is a 
nontrivial value depending on the state of the system. 

Now we first check the extensiveness of ^p N 2 \ The 
fit of the numerical data in Fig. 1 191 gives p( N ~ 2 ^ ~ jV - 99 . 

This tells us p( N ~ 2 ) (and consequently ^ p N ~ 2 ^) is pro- 
portional to the system size N, i.e., it is an extensive 
quantity. The dependence of p^ for the largest Lya- 
punov exponent on N is shown in the same plot. It is 
fitted with ~ TV - 76 , i.e., the Lyapunov vector for the 
largest LE is highly localized in space [Tfl I2H ■ This 
is consistent with our above observation (see Fig. [Bj. 

Then we study the temperature dependence of this 
newly defined length scale. For a 200-particle system, 
the temperature is increased from 0.4 to 1.6. Results of 
the simulation are shown in Fig. |2J3 where the variation 
of with temperature T is also presented. From the 
plot, one can see that p^ N ~^ (and consequently £p ) 
decreases gradually with the increase in temperature. 
This agrees with the intuitive expectation that, increas- 
ing of the temperature makes the fluctuation in the sys- 
tem stronger and stronger and renders the system be- 
coming more and more disordered. 

IV. CONCLUSION AND DISCUSSION 

In this paper, we presented numerical results about the 
Lyapunov instability of a Lennard-Jones system. Our 
simulations show that the step-wise structure found in 
the Lyapunov spectrum of hard-core systems disappears 
completely here. This is conjectured due to the strong 
fluctuations in the finite-time LEs |10| . A new technique 
|40| based on a spatial spectral analysis is employed to 
reveal the hidden long wave-length structure in LVs. A 
significantly sharp peak with low wave-number is found 
in the resultant spatial Fourier spectrum for LVs with 



A ~ 0. This serves as strong evidence that hydrodynamic 
Lyapunov modes do exist in the soft-potential systems 
|41j. Another important finding is that the dispersion 
relation for hydrodynamic Lyapunov modes, A^ Q ) versus 
kmax, appears to be universal irrespective of system tem- 
perature and particle density. 

In the study of two-dimensional and quasi-one- 
dimensional hard-core systems, two kinds of hydrody- 
namic Lyapunov modes are identified 0, 0|. One is 
referred to as transverse and the other is called longitudi- 
nal. The former modes do not propagate while the latter 
can [l^. According to this classification, for the trans- 
verse Lyapunov modes, taking the time average can be 
an useful way to identify the wave-like structure. In con- 
trast to this, the detection of the longitudinal Lyapunov 
modes is a relatively difficult task since due to its prop- 
agation time-averaging is no longer a suitable method to 
suppress the fluctuations . For the case of soft poten- 
tial system, strong fluctuations in local instabilities leads 
to the occasional mixing among Lyapunov vectors. This 
is partially reflected in the intermittent time evolution of 
spatial Fourier spectrum of LVs. Therefore the hydro- 
dynamic Lyapunov modes in a soft-potential system are 
more vague and more difficult to observe pl|. 

In the study of the hard-core case, it is conjectured 
that degeneracies in the Lyapunov spectrum and in wave- 
numbers of hydrodynamic Lyapunov modes are deter- 
mined by the intrinsic symmetries of the Hamiltonian and 
the boundary conditions. There is no contradiction be- 
tween this statement and the results reported here. The 
crucial point is the life-time of hydrodynamic Lyapunov 
modes. The above statement is for the ideal case of pure 
modes with infinite long life-time. For the Lennard-Jones 
system studied here, strong fluctuation in local instabil- 
ity leads to the mutual interaction and mixing among 
modes, which renders the life-time of modes becoming fi- 
nite. On the other hand, Lyapunov exponents are global 
quantities of the system, which are the result of a time- 
average along a long trajectory wandering all over the 
phase space allowed. Due to the mixing among modes, 
they no longer correspond to certain pure modes but to 
a mixture of several modes. Therefore the degeneration 
predicted on the basis of a symmetry analysis can not be 
seen here. Actually, fluctuations in local instabilities do 
exist in all dynamical systems. However, for the hard- 
core systems, it is relatively weak in the regime of A w 
(see Fig. 12 in [3 for the comparison of fluctuations in 
local instabilities for the two cases). There the mixing 
among modes is quite rare and weak. The life-times of 
the modes with A rj are so long that one can not feel 
it during the finite simulation time. Evidence for our 
argument above comes from the fact that as the fluctu- 
ations in local instabilities of hard-core system become 
stronger and stronger in increasing LEs from A « 0, the 
degeneration in LEs becomes worse and worse, i.e. the 
steps become steeper and steeper (see Fig. 12 and Fig. 8 
in Ref. 10,]). In this work, Fourier spectral analysis has 
been shown to be quite successful in detecting the hidden 
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wave-like structures in LVs. A more general theoretical 
consideration of this method will be given elsewhere |4jj . 

Until now, only the coordinate part of the LVs are used 
in the study of hydrodynamic Lyapunov modes. For the 
case of hard-core system, this is reasonable due to the 
interesting feature of LVs found in [12| that the angles 
between the coordinate part and the momentum part 
are always small, i.e, the two vectors are nearly paral- 
lel. For our soft potential system, we find that the angles 
between the coordinate part and the momentum part are 
no longer as small as in the hard-core systems. Even for 
the hard-core system, recent results show that the two 
vectors are not always parallel [2^| . Therefore it is neces- 
sary to reconsider the momentum part of LVs. Due to the 
specific feature of Hamiltonian systems, the momentum 
part of a LV with A > is identical to the coordinate part 
of a LV with —A. This ensures that the study of only the 
coordinate part of LVs is sufficient. But now one need to 
deal with LVs corresponding to both positive and nega- 
tive LEs. In general, the behavior of LVs with negative 
LEs should be different from those with positive ones. 
For the case reported in the current paper, there is no 
wave-like structure being observed in the former. 

We studied also the influence of density and temper- 
ature changes on the features of LVs and LEs. One ef- 
fect of decreasing the density is that the Lyapunov spec- 
trum becomes more and more bended. The relation of 



the bending in the Lyapunov spectrum with the sepa- 
ration of time scales was discussed recently in the work 
of Taniguchi and Morriss j2j|. It is obvious that, the 
collisions between particles become more and more rare 
as the density is decreased. They turn to be highly lo- 
calized events since they happen at only a few places at 
one moment. Therefore, the time scale of local collision 
events, which is related to the largest LEs |2|, is well sep- 
arated from that for the collective motion of the system, 
corresponding to the near-zero LEs. A further point is 
the sudden jump found in k max of the time-averaged spa- 
tial Fourier spectrum of LVs. It divides the whole set of 
LVs into two groups. The place of this sudden jump is 
in the regime where the Lyapunov spectrum is strongly 
bended. Further work is needed to reveal the underlying 
connection between these phenomena. 
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